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ABSTRACT 

We present the results from broadband X-ray spectral analysis of 3C 206 and PKS 0707—35 with Suzaku and 
Swift /BAT, two of the most luminous unobscured and obscured radio-loud active galactic nuclei (AGNs) with 
hard X-ray luminosities of 10 45 5 erg s 1 and 10 44 9 erg s~ 1 (14-195 keV), respectively. Based on the radio core 
luminosity, we estimate that the X-ray spectrum of 3C 206 contains a significant (~60% in the 14-195 keV band) 
contribution from the jet, while it is negligible in PKS 0707—35. We can successfully model the spectra with the jet 
component (for 3C 206), the transmitted emission, and two reflection components from the torus and the accretion 
disk. The reflection strengths from the torus are found to be R toms (=Q./2jr) = 0.29 ± 0.18 and 0.41 ± 0.18 for 3C 
206 and PKS 0707—35, respectively, which are smaller than those in typical Seyfert galaxies. Utilizing the torus 
model by Ikeda et al., we quantify the relation between the half-opening angle of a torus (0 oa ) and the equivalent 
width of an iron-K line. The observed equivalent width of 3C 206, <71 eV, constrains the column density in the 
equatorial plane to A^ q < 10 23 cm -2 , or the half-opening angle to 0 o a > 80° if A^ q = 10 24 cm~ 2 is assumed. That 
of PKS 0707—35, 72 ± 36 eV, is consistent with A^ q ~ 10 23 cm' 2 . Our results suggest that the tori in luminous 
radio-loud AGNs are only poorly developed. The trend is similar to that seen in radio-quiet AGNs, implying that 
the torus structure is not different between AGNs with jets and without jets. 
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Online-only material: color figures 


1. INTRODUCTION 

It has been recognized that feedback from active galactic 
nuclei (AGNs) plays a significant role in the coevolution 
between galaxies and supermassive black holes (SMBHs) in 
their centers (for a review see, e.g., Fabian 2012). AGN outflows 
in the form of jets and winds could expel and heat the gas in 
the interstellar medium, thus regulating both star formation and 
mass accretion onto SMBHs. The powerful jets from AGNs 
can affect the surrounding structure even at the scale of a 
galaxy cluster. Despite the importance of AGN feedback, several 
fundamental questions still remain open, such as: what are the 
physical mechanisms to launch relativistic jets and how the 
nuclear structure of AGNs with jets (radio-loud AGNs) differs 
from those without jets (radio-quiet AGNs)? 

X-ray observations bring us significant insight into the nuclear 

structure of AGNs. There are previous works that compare 
X-ray spectral properties of radio-loud and radio-quiet AGNs. 
Piconcelli et al. (2008) find rich emission lines in the soft 
X-ray spectrum of 3C 234, similar to those of radio-quiet AGNs. 
These lines come from a photoionized plasma, suggesting that 
the narrow-line region in radio-loud and radio-quiet AGNs has 
similar geometrical and physical properties. Warm absorbers 
similar to those found in radio-quite AGNs are also detected 
from several radio-loud AGNs (Torresi et al. 2012). Tombesi 
et al. (2010) discover ultra-fast outflows in some radio-loud 
AGNs, and their features are very similar to those previously 
observed from radio-quiet AGNs. However, no consensus has 
yet been firmly established whether or not there are differences 
between the two populations in their accretion disk structures 
and circumnuclear environments. To pursue these issues, our 


teams have been working in Suzaku follow-up of hard X-ray- 
selected radio-loud AGNs (Tazaki et al. 2010, 2011) detected 
in the Swift /BPsT survey (Cusumano et al. 2010; Baumgartner 
et al. 2013). 

The torus is a key structure in AGNs, which is considered 
to be responsible for supplying mass onto SMBHs. It causes 
obscuration observed in the so-called type-2 AGNs, a major 
AGN population in the universe (e.g., Ueda et al. 2003). To 
understand the origin of the torus, it is critical to identify 
the essential physical parameters that determine its structure. 
Previous studies, mainly based on the radio-quiet population, 
suggest that AGN luminosity is one of the most important 
parameters. The fact that the fraction of obscured AGN is anti- 
correlated with X-ray luminosity above Lx > 1 0 42 erg s' 1 (e.g., 
Ueda et al. 2003; Hasinger 2008; Burlon et al. 2011) indicates 
that the solid angle covered by the torus decreases with Lx- 
This correlation can also explain the so-called X-ray Baldwin 
effect (Iwasawa & Taniguchi 1993) that the equivalent width of 
the iron-Ka line decreases with X-ray luminosity (Grandi et al. 
2006; Ricci et al. 2013). 

In this paper, we study the spectra of two very luminous 
radio-loud AGNs, 3C 206 and PKS 0707-35, observed with 
Suzaku and Swf/f/BAT to reveal their nuclear structure. The 
main purpose is to establish if the correlation of the torus 
opening angle with luminosity seen in radio-quite AGNs also 
holds for radio-loud AGNs. For quantitative discussion, we 
utilize the numerical torus model by Ikeda et al. (2009) to 
predict the equivalent width of the iron-K line as a function 
of torus opening angle and column density. 

Our targets, 3C 206 and PKS 0707—35, are the most luminous 
unobscured (type-1) and obscured (type-2) radio-loud AGNs 
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Table 1 

List of Target 


Tazaki et al. 


Target Name ( Swift ID) 

3C 206 (J0839.6 1213) 

PKS 0707-35 (J0709.4— 3559) 

R.A. (J2000) a 

08 39 50.62 

07 09 14.09 

Deck (J2000) a 

-12 14 33.9 

-36 01 21.8 

Redshift 

0.198 

0.111 

Suzaku observation ID 

705007010 

706008010 

Nominal position 

HXD 

XIS 

Start time (UT) 

2010-05-08T16:32:34 

201 l-04-02T02:16:10 

End time (UT) 

2010-05- 10T1 1:45: 10 

201 l-04-04T01:32:13 

Exposure 8 (XIS) (ks) 

80.6 

81.1 

Exposure (HXD/PIN) (ks) 

65.6 

57.2 


Notes. 

a The position of each source is taken from the NASA/IPAC Extragalactic Database. 
b Based on the good time interval of XIS-0. 


having hard X-ray luminosities of £ 14—195 lev = 10 455 erg s _1 
and 10 44 ' 9 erg s _1 , respectively, among non-blazar type AGNs 
with fluxes above 2 x I 0“ 1 1 erg crrT 2 s ” 1 (14-195 keV) detected 
in the Swift/ BAT 58 month survey (Baumgartner et al. 2013). 
3C 206 (z = 0.1979; Ho & Kim 2009) is a radio-loud quasar, 
classified as Seyfert 1.2 from the optical spectrum (Veron- 
Cetty & Veron 2006). The source is located in a cluster of 
galaxies with the Abell richness class 1 (Ellingson et al. 1989). 
The observation with the EINSTEIN High Resolution Imager 
detected no significant extended X-ray emission within ~2 
arcmin from the AGN (Hall et al. 1997). The host galaxy of 
3C 206 has an elliptical morphology. The black hole mass is 
estimated to be M BH = 10 8 - 8-8 - 9 M e (Sikora et al. 2007; Liu 
et al. 2006) from the empirical formula using the continuum 
luminosity at 5100 A and the line widths of Ha and H/i. PKS 
0707—35 (z = 0.1108; Burgess & Hunstead 2006) is a radio 
galaxy of Fanaroff-Riley type 11. This object has a largely 
extended radio structure up to ~500". The galaxy morphology 
and black hole mass of PKS 0707—35 are not known. We 
classify this source as an obscured (type- 2 ) object on the basis 
of the X-ray spectrum, which shows an intrinsic absorption with 
a hydrogen column density of Ah > 10 22 cm 2 (Section 3.2.2). 

This paper is organized as follows. Section 2 summarizes the 
observations. In Section 3, we present the data analysis and the 
results including the light curves of Suzaku and detailed model 
fit to the broadband spectra obtained with Suzaku and Swift/ 
BAT. A summary and discussion are given in Section 4. The 
cosmological parameters ( Hq , Q m , il,_) — (70 km s -1 Mpc -1 , 
0.3, 0.7; Komatsu et al. 2009) are adopted in calculating the 
luminosities. The errors on the spectral parameters correspond 
to the 90% confidence limits for a single parameter. 

2. OBSERVATIONS 

We observed the two radio-loud AGNs 3C 206 and PKS 
0707-35 with Suzaku on 2010 May 8-10 and 201 1 April 2-4, 
respectively, for a net exposure of ~80 ks. The hard X-ray 
fluxes in the 14-195 keV band averaged from 2005 to 2010 
(Baumgartner et al. 20 1 3) 5 are 2.5 x 10~ u erg cm ” 2 s -1 (3C 
206) and 2.2 x KT 11 erg cm ” 2 s ” 1 (PKS 0707-35), which 
are sufficiently bright for detailed broadband spectral analysis. 
The basic information of these targets and observation logs are 
summarized in Table 1 . 

Suzaku, the fifth Japanese X-ray satellite (Mitsuda et al. 
2007), carries three currently active X-ray CCD cameras called 

5 http://heasarc.gsfc.nasa.gov/docs/swift/results/bs58mon/ 


the X-ray Imaging Spectrometer (two Front-side Illuminated 
ones (XIS-FIs), XIS-0, 3, and one Back-side Illuminated one 
(BI-XIS), XIS- 1 ) and the Hard X-ray Detector (HXD) composed 
of Si PIN photodiodes and gadolinium silicon oxide (GSO) 
scintillation counters. We utilize the data of the XIS-FIs, XIS-BI, 
and HXD /PIN in the energy bands of 1-12 keV, 0.5-8 keV, and 
16-35 keV (PKS 0707—35 only), respectively, where the source 
detection is not affected by the background uncertainties and the 
calibration is reliable. We do not use the HXD/PIN data of 3C 
206 because of the unusually high count rate of the non-X-ray 
background (NXB) that prevented a firm detection of the source 
signal. 

The data of HXD/GSO, covering above 50 keV, are not 
utilized, because the fluxes of both targets are too faint to 
be detected. 3C 206 and PKS 0707—35 were focused on the 
nominal center position of the HXD and the XIS, respectively. 
For spectral analysis, we also include the Swift / BAT spectra 
covering the 14-195 keV band averaged over 58 months from 
2005 to 2010 (Baumgartner et al. 2013). 

3. ANALYSIS AND RESULTS 

We reduce the unfiltered event files processed with the ver- 
sions of 2.5.16.28 (3C 206) and 2.5.16.29 (PKS 0707-35). To 
reflect the latest energy calibration, we apply the xispi ftool with 
the CALDB files ae_xi{0,l,3}_makepi_201 10621. fits. These 
events files are then screened with standard selection criteria 
described in the Suzaku ABC Guide. To extract light curves 
and spectra, we use the FTOOLS package (heasoft version 
6.11). Spectral fitting is performed with the XSPEC software 
(version 12.7.0o). The XIS events of each target are extracted 
from a circular region centered on the source peak with a ra- 
dius of 3', and the background is taken from a source-free re- 
gion within the XIS field of view. We generate the redistri- 
bution matrix file (RMF) and ancillary response file (ARF) of 
the XIS with the xisrmfgen and xissimarfgen ftools (Ishisaki 
et al. 2007), respectively. The “tuned” NXB event files pro- 
vided by the HXD team are utilized for the background sub- 
traction of the HXD/PIN data. The simulated cosmic X-ray 
background is added to the NXB spectrum based on the for- 
mula given by Gruber et al. (1999). We use the HXD/PIN re- 
sponse file ae_hxd_pinxinome9_20100731.rsp. The spectra of 
Swift / BAT and response matrix are available from the 
HEASARC Web site of the 58 month survey . 6 


6 http://heasarc.gsfc.nasa.gov/docs/swift/results/bs58mon/ 
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Figure 1. (a) Light curves of 3C 206 in the 0.5-2 keV (upper) and 2-10 keV 
bands (middle). The lower panel (RATIO 1) shows the hardness ratio of 
2-10 keV/0.5-2 keV. (b) Light curves of PKS 0707-35 in the 0.5-2 keV, 
2-10 keV, and 16-35 keV bands. The lower two panels (RATIO 1 and RATI02) 
show the hardness ratios of 2-10 keV/0.5-2 keV and 16-35 keV/2-10 keV, 
respectively. 


3C 206 



PKS 0707-35 



(b) 

Figure 2. (a) Folded spectra of 3C 206 observed with the Suzaku XIS-FIs (black 
crosses), XIS-BI (red crosses with open triangles), and Swiyi/BAT (blue crosses 
with open squares). The normalizations of the XIS-BI spectrum are shifted for 
clarity of the plots, (b) The same for PKS 0707—35. The Suzaku PIN spectrum 
is included (magenta crosses with open circles). 

(A color version of this figure is available in the online journal.) 


3.1. Light Curves 

Figure 1(a) shows the light curve of 3C 206 in the 0.5-2 keV 
band obtained from X1S-1 (upper), that in the 2-10 keV band 
combined from XIS-0 and XIS-3 (middle), and the hardness 
ratio between the two bands (RATIO 1). Figure 1(b) represents 
those of PKS 0707—35, together with the HXD/PIN light curve 
in the 16-35 keV band. The hardness ratio between the 2-10 keV 
and 16-35 keV bands (RATI02) is also plotted. The bin width 
is 5760 s, the orbital period of Suzaku, chosen to eliminate any 
systematics caused by the orbital change of the satellite. The 
background is subtracted. As noted from the figures, we find 
no significant time variability in RATIO 1 and RATI02. Hence, 
we analyze the time-averaged spectra over the whole Suzaku 
observations for both objects. 

3.2. Spectral Analysis 

In the spectral analysis, we simultaneously fit the Suzaku and 
Swift /BAT spectra. The XIS data in the 1.7-1. 9 keV band are 
excluded because of calibration uncertainties associated with 


the instrumental Si-K edge. Figure 2 shows the spectra of 3C 
206 and PKS 0707—35 folded with the detector response. The 
relative flux normalization of the XIS-BI to the XIS-FI is set to 
be a free parameter, while that of the HXD/PIN to the XIS-FIs 
is fixed at 1.16 for PKS 0707—35 (XIS nominal position), based 
on the calibration using the Crab Nebula (Maeda et al. 2008). 
The Galactic absorption, /V^' 11 ' , is always included in the spectral 
model with phabs in the XSPEC terminology, which is fixed 
at the value derived from the Hi map of Kalberla et al. (2005; 
Ag al = 5.68 x 10 20 cm" 2 for 3C 206 and 1.70 x 10 21 cm’ 2 
for PKS 0707—35). Solar abundances by Anders & Grevesse 
(1989) are assumed throughout our analysis. 

We model the emission from the nucleus (other than the 
jet component) with a power law with exponential cutoff 
(E 1 x exp(— E/E cal )) and Compton reflection components 
from cold matter. Unlike most of previous works, in this paper 
we consider two reflection components, one from the accretion 
disk and the other from the torus, in order to achieve a physically 
consistent picture. Since the direct emission is subject to time 
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variability between the Suzaku and Swift /BAT observations, we 
make the relative normalization of the direct component and its 
disk reflection (NormBAi) free between the two spectra. Since 
we do not anticipate significant time variability in the reflec- 
tion component from the torus, which usually has a size of an 
order of ~0.1 pc or larger (Suganuma et al. 2006), the rela- 
tive normalization of the torus reflection is tied between the 
Suzaku and Swift/ BAT observations. In fact, when we apply 
the formula by Ghisellini & Tavecchio (2008), the torus dis- 
tances of 3C 206 and PKS 0707—35 are estimated to be ~2 pc by 
assuming a bolometric correction factor of Lb 0 i/T 2 — tokev = 30 
for the disk luminosity. For both targets, the cutoff energy is 
fixed at 200 keV, which is suggested from several radio-loud 
AGNs (e.g., Grandi et al. 2006; Sambruna et al. 2009). The av- 
eraged cutoff energy for local Seyfert galaxies is reported to be 
~300 keV by Dadina (2008), which is somewhat higher than 
that of radio-loud AGNs. We confirm that the results are little 
affected even if E c ut = 300 keV is adopted instead. 

To calculate the reflection components, we adopt the pexmon 
model (Nandra et al. 2007) in XSPEC, where the fluorescence 
lines of iron-Ka, iron-K/J, and nickel-Ka are self-consistently 
included in the reflected continuum computed by the pexrav 
code (Magdziarz & Zdziarski 1995). pexrav model assumes 
a neutral, Compton-thick (Ah > 10 24 cm -2 ) medium with a 
plane geometry. The relative strength to the direct component 
is defined by R = Q./2jt, where Q is the solid angle of the 
reflector viewed from the central irradiating source. We consider 
the Doppler and relativistic smearing in the disk-reflection 
component, which is thus modeled as rdblur * pexmon in 
the XSPEC terminology. Since it is difficult to constrain the 
parameters of rdblur from our spectra, we assume the inner 
and outer radii of the accretion disk to be 1 00 r g (r g = GM/c 2 
is the gravitational radius) and 10 5 r g , respectively, with a typical 
emissivity law of r _3 . The adopted inner radius can be regarded 
as a typical value observed from radio galaxies (Tazaki et al. 
2010; Sambruna et al. 2009; Larsson et al. 2008), although it is 
yet uncertain. The inclination angles between the normal axis 
of the accretion disk and the line of sight are assumed to be 
30° for 3C 206 (type-1) and 60° for PKS 0707-35 (type-2). 
Even when we adopt an inclination angle of 18° for 3C 206 
(the lowest available value in the pexrav/pexmon model), the 
result of broadband fitting is essentially unaffected. Such a 
low inclination is unlikely in our case, however, because the 
predicted jet emission largely overpredicts the observed X-ray 
flux (see Section 3.2.1 for the details). As for the reflection 
component from the torus, we ignore any smearing effects, 
which cannot be resolved at the CCD energy resolution. It is 
described as a single pexmon model, whose inclination angle 
is fixed at 60° for both targets, as a representable value for the 
complicated torus geometry. 

3.2.1. Broadband Spectral Analysis of 3C 206 

Fit without jet component. To evaluate possible contamination 
of X-ray emission from the galaxy cluster around the AGN, we 
compare the radial profile of 3C 206 obtained from the XIS-0 
and XIS-1 images in the 0.5-2 keV band with that of Mrk 205, an 
AGN not located in a cluster, observed in 2010 May at the HXD 
nominal position ( Suzaku observation ID: 705062010). The 
results are plotted in Figure 3. There is no evidence for diffuse 
emission from the cluster of galaxies surrounding 3C 206, which 
could be extended over ~600 kpc (~3' at z — 0.1979). This is 
consistent with the result by Hall et al. (1997) using Einstein. We 



0 12 3 


arcmin 

Figure 3. Background-subtracted radial profiles of 3C 206 (black crosses with 
open circles) and Mrk 205 (red crosses with open squares) obtained from the 
XIS-0 (upper) and XIS-1 (lower) images in the 0.5-2 keV band. 

(A color version of this figure is available in the online journal.) 

thus analyze the spectra of 3C 206 by neglecting any emission 
from the cluster of galaxies. 

For the spectral fit of 3C 206, we utilize the data of XIS- 
FIs in the 1-12 keV band, XIS-BI in the 1-8 keV band, and 
Swift/ BAT in the 14-195 keV band; because the XIS spectra 
show a soft excess below 1 ke V over a power law extrapolated 
from higher energies, we do not include the XIS-BI spectrum 
below 1 keV in the broadband spectral analysis to avoid com- 
plexity. Possible origins of the soft excess will be discussed later 
in this subsubsection. First, we apply a simple model consisting 
of a direct continuum and its reflection component only from the 
torus, expressed as “phabs (constant * zpowerlw * zhighect 1 + 
pexmon) ” in the XSPEC terminology. Here, the “constant” fac- 
tor is introduced to represent the time variability between the 
Suzaku and Burst Alert Telescope (BAT) observations. We ob- 
tain the best-fit photon index of T = 1.75 ± 0.01 and reflec- 
tion strength from the torus of R to ms = 0.25 ± 0.09 relative to 
the time-averaged flux of the direct component determined by 
Swift /EKE with x 2 /dof = 829.7/825. 

Next, we adopt a model expressed as “phabs (constant 
(zpowerlw * zhighect + rdblur * pexmon) + pexmon)”, 
adding the reflection component from the accretion disk to the 
previous model to achieve a physically consistent picture, even 
though the previous fit without the disk reflection is acceptable 
in terms of the y 1 value. Since it is difficult to determine the 
strengths of the two reflection components independently from 
our data, we fix only that of the disk component. According 
to the calculation by George & Fabian (1991), the equivalent 
width of the iron-K line is predicted to be EWdisk = 17-31 eV 
for an accretion disk truncated at 1 00 r„ for a power-law 
photon index of 1.7-1. 9. Here, we assume an inclination 
angle of 30° and scale height of the irradiating source of 10 
We then estimate the corresponding reflection strength of 
R(= Q/2tt) ~ 0.1 in the pexmon model so that the predicted 
equivalent width (17-31 eV) is consistently reproduced. Thus, 
we fix R = 0.1 for the reflection component from the accretion 
disk. We find that the Suzaku and BAT spectra of 3C 206 


7 In the zhighect model, we set the threshold energy above which the 
exponential cutoff is applied to be zero. 
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favor this model with /?dj s k = 0.1 (x 2 /dof = 825.3/825) 
over the single reflection model, although the difference in the 
X 1 value, A/ 2 = 4.4, is not highly significant. We obtain 
the best-fit photon index of T = 1.75±0.01 and reflection 
strength from the torus of R toms = 0.16 ± 0.09 relative to 
the time-averaged flux of the direct component determined by 
Swift/ BAT. 

Fit with jet component. Emission from the jets may be 
contained in the X-ray spectra of radio-loud AGNs. To estimate 
the jet contribution in the X-ray band, we utilize the radio 
luminosity in the AGN core, which is attributable to the emission 
from the innermost region of the jet. According to the unified 
scheme of AGNs (Antonucci 1993), radio-loud quasars and 
radio galaxies are intrinsically the same populations of blazars 
except that the jet axis is not aligned with the line of sight. 
Thus, from the observed radio luminosity, we can estimate 
the broad band spectral energy distribution (SED) of the jet 
emission in (non-blazar) radio-loud AGNs, by referring to 
averaged SEDs of blazars given as a function of radio luminosity 
(the so-called blazar sequence) after correcting for relativistic 
beaming. The observed SED from an emitter moving with a 
relativistic velocity of v shifts toward higher frequencies by 
a beaming factor of S (= l/y(l — f cos0), where f = v/c, 
y = \/yJ\ - f 1 , and 0 is the angular separation from the line 
of sight), and toward higher luminosities by <5 4 (Section 5.8. 1.1 
in Beckmann & Shrader 2012). By assuming the inclination 
angle of Q mc = 30° and the jet velocity of j> = 0.9-0.99 
for type-1 radio-loud AGNs, the beaming factor of 3C 206 
is estimated to be /> 3 c 206 = 0.5-2. Typically, blazar ~ 10 in 
blazars (Dondi & Ghisellini 1995). Then, we can convert the 
SED of typical blazars into that of 3C 206 by multiplying a 
factor of <5 3C2 06/<5biazai- and (<$ 3C 206/<5biazar) 4 in the frequency and 
luminosity, respectively. 

Using the blazar SEDs given in Donato et al. (2001) and 
the radio luminosity of 3C 206 (vL ,, = 2.3 x 10 42 erg s -1 at 
5 GHz; Elvis et al. 1994), we find that the intrinsic SED of the 
jet emission in 3C 206 corresponds to those of the three most 
luminous blazar classes in Figure 10 of Donato et al. (2001). 
We estimate that the observed X-ray luminosity from the jet in 
3C 206 should be vL v ~ io 4L7 ~ 42 - 8 erg s” 1 at 1 keV, which 
is ~l%-20% of the total observed flux, and the spectrum is 
approximated by a power law with a photon index of 1.2- 1.4 
in the X-ray band below 100 keV. Note that when we assume 
an inclination angle of 9 mc = 18° (<5 3 c206 = 1-7-2), the jet 
luminosity is estimated to be vL v ~ 1 0 45 1 44 2 erg s _1 at 1 keV, 
which is >10 times higher than the total observed flux. 

To take into account the jet emission, we add a power law 
to the spectral model described above. Considering a large 
uncertainty in the predicted luminosity of the jet component, 
however, we fix the photon index at 1 .3 but set the normalization 
to be a free parameter. For simplicity, we also assume that the 
flux ratio between the jet (power law) and the direct component 
from the disk (cutoff power law) did not change between 
the Suzciku and Swift /HAT observations. The final model, 
which is shown in Figure 4, is expressed as “phabs (constant 
(zpowerlw * zhighect + rdblur * pexmon + zpowerlw) + 
pexmon)”, where the constant factor takes account of flux 
variation between the two observations, and the first to fourth 
terms represent (1) the direct component, (2) the reflection from 
the accretion disk, (3) the jet component, and (4) that from the 
torus, respectively. 

We confirm that this model, including the jet emission, well 
reproduces the broadband spectrum in the 1-200 keV band 


3C 206 



Energy (keV) 


Figure 4. Unfolded spectra of 3C 206 obtained with the XIS-FIs in the 1-12 keV 
band (black crosses without mark) and Swift/ BAT in 14-195 keV band (blue 
crosses with open squares) in units of EF(E) ( F(E ) is the energy flux at the 
energy E). The XIS-BI spectrum is not plotted here but is included in the fitting. 
The top solid curves show the best-fit model (Table 2). The dotted (black), thin 
solid (red), thick solid (blue), and triple-dot-dashed (orange) curves correspond 
to the direct component, reflection component from the torus, that from the 
accretion disk, and jet component, respectively. The residuals in units of cr are 
plotted in the lower panel. 

(A color version of this figure is available in the online journal.) 

(X 2 /dof = 818.0/824). Besides the requirement from the con- 
sistency with the radio core luminosity, we confirm that the im- 
provement of the fit by adding the jet component is marginally 
significant with an T-test probability of 99.3%. It is found that 
the fraction of the jet flux in the total one is ~20% in the 
2-10 keV band and ~60% in the 14-195 keV band, respec- 
tively. This result is consistent with our rough estimate using 
the radio luminosity and blazar sequence. As for the non-jet 
emission, we obtain the best-fit photon index of T = 1.86 4 qq 8 7 
and reflection strength from the torus of R t orus = 0.29 ± 0.18. 
The best-fit parameters are summarized in Table 2. Although the 
improvement from the previous fit without the jet component 
is only marginally significant (at 99.3% confidence level from 
the T-test), we regard this model as the most realistic descrip- 
tion of the broadband data of 3C 206, considering that the jet 
contribution in the X-ray spectrum is quite likely. 

We check possible systematic errors caused by our assump- 
tion on the reflection component from the accretion disk, which 
is quite uncertain. If we assume 10 r g as the inner radius of 
the accretion disk instead of r m = 100 r g , we expect a stronger 
reflection of R^w ~ 0.5, which predicts EWdisk = 100-130 eV, 
based on the calculation by George & Fabian (1991). Adopting 
r ln = 10 r g and R^sk = 0.5 in the above model (including the 
jet component) yields the best-fit photon index of T = 1 .88/y q 7 6 
and reflection strength from the torus of R toms < 0.27 with 
X 2 /dof = 822.7/824. The fit becomes considerably worse but 
these parameters are still consistent with those obtained in the 
case of r m — 1 00 r„ and /Gist =0.1 within the statistical errors 
at the 90% confidence level. 

Soft excess. As mentioned above, we detect the soft excess 
below 1 keV in the XIS-BI spectrum when the above model is 
extrapolated toward lower energies; when the XIS-BI spectrum 
in the 0.6-1 keV band is included, the goodness of the fit 
becomes rather poor with x 2 /dof = 989.6/926. Similar features 
are observed in the X-ray spectra of many Seyfert 1 galaxies. 
First, we consider possible contribution of optically thin thermal 
emission from a gas unrelated to the AGN. By adding the 
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Table 2 

Best-fit Parameters of Broadband Spectra 




3C 206 

PKS 0707-35 

( 1 ) 

/Vg 31 ( 10 22 crrr 2 ) 

0.0568 a 

0.1 70 a 

( 2 ) 

Normal 

1.02 ± 0.01 

0.97 ± 0.02 

(3) 

NormpiN 


1.16 a 

(4) 

NomiBAT 

0 73 + 0.17 
U ‘ ' — 0.14 

0.833 1 * 

(5) 

A nuc- (photon cm -2 s -1 keV -1 ) 

3.40^o 2 i8 x 10 ~ 3 

1-79/q2o x 10~ 3 

( 6 ) 

Aj et (photon cm -2 s -1 keV -1 ) 

4.21/j 54 x 10-4 


(V) 

N £ ( 10 22 cm" 2 ) 


5 1 +3-2 
J -2.3 

( 8 ) 

IV 2 (10 22 cm" 2 ) 


2 O +0 - 4 
z ‘ u -0.9 

(9) 

/ 


n 41 + 0.43 
u>fl - 0.21 

( 10 ) 

r 

1 o/c+0.08 

1 ° -0.07 

1 66 +007 
1-0 —0.06 

(ID 

Gut (keV) 

200 a 

200 a 

( 12 ) 

^disk 

( 0 . 1 a ) 

( 0 . 2 a ) 

(13) 

n n ( r g) 

100 a 

100 a 

(14) 

^torus 

0.29 ± 0.18 (0.21/° 0 12 0 ) 

0.41 ±0.18 (0.34 ±0.12) 

(15) 

/scat (%) 


1.3 ±1.2 (1.0 ±0.9) 

(16) 

Tjet 

1.3 a 


(17) 

G -10 (erg cm -2 s -1 ) 

1.0 x 10 ~ n ( 2.6 x 10 ” 12 ) 

5.6 x 10 ~ 12 

(18) 

Go — 50 (erg cm ' 2 s -1 ) 

1.8 x 10 ~ n ( 8.0 x 10 ~ 12 ) 

1.6 x 10 ~ n 

(19) 

Z. 2-10 (erg s _1 ) 

2.8 x 10 44 (6.5 x 10 43 ) 

2.1 x 10 44 


X 2 /dof 

818.02/824 

616.31/582 


Notes. The errors correspond to 90% confidence level for a single parameter. Both/ sca t and R are normalized to 
the flux of the direct component as measured with Swift/BXT ( Suzaku ). 

(1) The hydrogen column density of Galactic absorption by Kalberla et al. (2005). 

(2) Normalization ratio between the XIS-BI and XIS-FI spectra. 

(3) Normalization ratio between the HXD/PIN and XIS-FI spectra. 

(4) Normalization ratio of the direct component between the BAT and Suzaku spectra. 

(5) The power-law normalization of the nuclear emission at 1 keV. 

(6) The power-law normalization of the jet emission at 1 keV. 

(7) Hydrogen column density for the direct component, which covers (100 x /)% in the line of sight. 

(8) Hydrogen column density for the direct component, which covers (100 x (1 — /))% in the line of sight. 

(9) Covering fraction. 

(10) The power-law photon index of the direct component. 

(11) The cutoff energy. 

(12) The strength of the reflection component from the accretion disk relative to the flux of the direct component 
as measured with Swift /BAT (Suzaku), defined as R = Q./2n, where Q is the solid angle of the reflector. 

(13) The assumed inner radius of the accretion disk. 

(14) The strength of the reflection component from the torus relative to the flux of the direct component as 
measured with Swift /BAT (Suzaku), defined as R = Q./2n, where Q. is the solid angle of the reflector. 

(15) The scattered fraction relative to the flux of the direct component as measured with Swift /BA lT (Suzaku). 

(16) The power-law photon index of the jet component. 

(17) The observed Suzaku flux (that of the jet component only) in the 2-10 keV band. 

(18) The observed Suzaku flux (that of the jet component only) in the 10-50 keV band. 

(19) The 2-10 keV intrinsic luminosity (that of the jet component only) obtained with Suzaku corrected for the 
absorption. 

a The parameters are fixed. 


apec model in the XSPEC terminology, the fit is significantly 
improved (/ 2 /dof = 911.1/924) with a temperature of kT ~ 
0.15 keV and a luminosity of ~ 4 x 10 43 erg s _1 . 

However, this luminosity is too high as that from star-forming 
regions in the host galaxy of AGNs (e.g., Cappi et al. 1999). As 
mentioned in the previous section, this could not be the hot gas 
in the galaxy cluster, either. Thus, we conclude that this model 
is unrealistic, even though it is statistically acceptable. 

A possible origin is a thermal Comptonization component 
of soft seed photons from the disk (e.g., Noda et al. 2011), 
although the physical picture is not yet established. We find that 
the compTT model (Titarchuk 1994) successfully represents 
the soft excess (x 2 /dof = 902.6/923 corresponding to an 
/'’-test probability of > 99.99%) with the best-fit electron 
temperature of kT e ~ 42 keV and optical depth of r ~ 0.1 for 


blackbody seed photons of 0.01 keV. The temperature of this 
Comptonization corona is consistent with the result from the 
Seyfert 1 galaxy MCG 6-30-15 obtained by Noda et al. (201 1), 
while the optical depth is smaller. We need to be aware of the 
degeneracy between the temperature and optical depth, however. 
The observed flux of this component is 5 x 10~ 13 erg cm~ 2 s' 1 
in the 0.5-1 keV band. 

An alternative explanation of the soft excess is reflection 
from a mildly ionized disk, which emits strong iron-L lines 
below 1 keV (e.g., Ross & Fabian 1993; Crummy et al. 2006; 
Sambruna et al. 2011). Accordingly, we utilize the reflionx 
model (Ross & Fabian 2005; Ross et al. 1999) convolved with 
rdblur (Fabian et al. 1989), instead of “cold” reflection modeled 
as rdblur * pexmon. Since such an ionized disk region is 
expected to be located close to the SMBH, we fix r m = 6 r g , the 
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innermost stable circular orbit of a non-rotating black hole. This 
model also reproduces the soft excess of our spectra with the 
best-fit ionization parameter / ~ 200 erg cm -1 s _1 , which 
predicts prominent iron-L line features below 1 keV (e.g., 
Ross & Fabian 1993). We obtain x 2 /dof = 91 1.0/923 with an 
/•’-test probability > 99.99% from the model without the soft- 
excess component. The reflection strength roughly corresponds 
to R ~ 0.4 as estimated from its continuum flux above 10 keV. 
While we cannot distinguish which is the more realistic model 
explaining the soft excess, our results on the reflection from the 
torus are not affected. 

3.2.2. Broadband Spectral Analysis ofPKS 0707—35 

In the spectral analysis of PKS 0707—35, we use the data of 
the XIS-FIs in the 1.0-12 keV band, XIS-BI in the 0.6-8 keV 
band, HXD/PIN in the 16-35 keV band, and Swift /BAT in the 
14-100 keV band. Swift /BAT did not detect positive signals 
above 100 keV from this source. As described below, the jet 
component is negligible in the X-ray spectrum ofPKS 0707—35. 
We find the continuum to be subject to an intrinsic absorption 
with Nh ~ 3 x 10 22 cm~ 2 , on which we base our classification of 
this source as an obscured AGN. Previous studies (e.g., Turner 
et al. 1997) showed that a scattered (or photoionized) component 
of the direct emission by a gas surrounding the nucleus is 
almost always observed in the X-ray spectra of absorbed AGNs, 
including radio-loud ones (e.g., Piconcelli et al. 2008; Tazaki 
et al. 2011; Grandi et al. 2007; Torresi et al. 2009; Evans et al. 
2010; Balmaverde et al. 2012). We thus essentially adopt the 
same model components as applied to 3C 206 except that (1) 
no jet component is included, (2) intrinsic neutral absorption is 
applied to the model components, and (3) a scattered component 
is added with the same photon index as that of the transmitted 
one. Regarding the third component, we do not find significant 
improvement of the fit by adding emission lines expected from 
a photoionized plasma to the power-law model due to limited 
statistics and energy resolution of the data. 

We confirm the absence of a significant jet flux in the X-ray 
spectrum of PKS 0707—35 based on the same discussion as for 
3C 206. The core radio flux of PKS 0707—35 is measured to 
be 25 mJy at 2.3 GHz by lones et al. (1994). Assuming the 
inclination angle of 0 mc = 60°, which is larger than the case 
of 3C 206 because it is a type-2 AGN, and the jet velocity of 
P — 0.9-0.99, we estimate the beaming factor of the jet in PKS 
0707—35 to be <5 = 0.1-1. By comparison with the SEDs of 
blazars (Donato et al. 2001), we estimate that the fraction of the 
jet flux is less than ~0.3% of the total flux at 1 keV corrected 
for the absorption. 

As in the spectral analysis of 3C 206, we first consider the 
reflection component only from the toms, adopting a model ex- 
pressed as “phabs (constant * zphabs * zpowerlw * zhighect + 
zphabs * pexmon + / scat * zpowerlw * zhighect)” in the 
XSPEC terminology. These terms correspond to the direct com- 
ponent, the reflection component from the torus, and the scat- 
tered component, respectively. The constant factor represents 
the flux variation between the Suzaku and BAT observations. 
We assume that the fluxes of the scattered and torus-reflection 
components did not vary between the two observations. Since it 
is difficult to constrain the intrinsic absorption separately for the 
direct and reflection components that are emitted from different 
regions, we adopt a same Ah value for them. We find that the 
model can reproduce the broadband spectra of PKS 0707—35 
reasonably well (x 2 /dof = 631.0/584) with T = 1.55 ± 0.04, 
Ah = (2.5 ± 0.1) x 10 22 cm -2 , R to rus = 0.61 ± 0.22, and 


PKS 0707-35 



Energy (keV) 


Figure 5. Unfolded spectra of PKS 0707—35 obtained with the XIS-FIs in the 
1-12 keV band (black crosses without mark), the HXD/PIN in the 16-35 keV 
band (magenta crosses with open circles), and Swift / BAT in 14-100 keV band 
(blue crosses with open squares) in units of EF(E). The top solid curves show 
the best-fit model (Table 2). The dotted (black), thin solid (red), thick solid 
(blue), and dot-dashed (magenta) curves correspond to the direct component, 
reflection component from the torus, that from the accretion disk, and scattered 
component, respectively. The residuals in units of o are plotted in the lower 
panel. 

(A color version of this figure is available in the online journal.) 

/ scat = (4 ± 1)% (the latter two parameters are defined as those 
relative to the time-averaged flux of the direct component as 
measured with Swift /BAT). 

Next, we add a reflection component from the accretion disk. 
The model is thus expressed as “phabs (constant * zphabs 
(zpowerlw * zhighect + rdblur * pexmon) + zphabs * 
pexmon + / scat * zpowerlw * zhighect)” in the XSPEC 
terminology, where the second term corresponds to the disk- 
reflection component. We estimate the strength of the reflection 
component from the accretion disk by assuming the inner radius 
of 100 and the inclination angle of 60° (instead of 30° as 
assumed for 3C 206). For a photon index of 1.5-1. 8, we find 
Rdisk ~ 0.2 and EWd; s k = 13-26 eV from George & Fabian 
(1991), and hence we fix Rdisk — 0-2. We find that the fit is 
better (x 2 /dof = 629.2/584) than the previous one with the same 
degree of freedom, although the significance of the improvement 
(Ax 2 = 1 . 8) is weaker than the case of 3C 206. As expected, we 
obtain a weaker reflection strength from the torus, R tom s ~ 0.5, 
while the other parameters are little changed. 

Since the above fit is not yet statistically acceptable, we then 
adopt a partial covering model instead of a single absorption, 
where two different absorbing column densities are considered 
with fractions / and 1 — /. This model often better reproduces 
the X-ray spectra of type-2 AGNs (see, e.g., Eguchi et al. 2009). 
The model, which is shown in Figure 5, is thus expressed as 
“phabs (constant * zphabs * zpcfabs (zpowerlw * zhighect + 
rdblur * pexmon) + zphabs * zpcfabs * pexmon + 
/scat * zpowerlw * zhighect)”. We find that the fit is sig- 
nificantly improved (x 2 /dof = 616.3/582) with an /-’-test prob- 
ability of 99.76%. The parameters of the partial covering are 
determined as A^ = 5.1/ 3 2 ' 2 3 x 10 22 cm~ 2 (41^%) and Ay = 
2.0/°q 4 9 x 10 22 cm -2 (59/42%). We obtain the best-fit photon 
index of T = 1 -66/006’ the reflection strength from the torus 
of R toms = 0.41 ±0.18 (to the BAT flux), and the scattering 
fraction of / scat = (1.3 ± 1.2)% (to the BAT flux). The best-fit 
parameters are summarized in Table 2. We regard this model as 
the best description of the spectrum of PKS 0707—35. 
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Finally, we also check systematic uncertainties related to 
the assumption of the reflection component from the accretion 
disk. If we assume the inner radius of the accretion disk to 
be 10 /‘ g , the calculation by George & Fabian (1991) predicts 
7? disk ~ 0.9 and EWd; s k = 80-110 eV for an inclination of 
60° and a photon index of 1.5-1. 8. Adopting /Gist = 0.9 and 
r; n = 1 0 r g , we obtain the following best-fit parameters with 
X 2 /dof = 614.6/582: T = 1.79 ± 0.06, R toms = 0.44 ± 0.19 
(to the BAT flux), / sca , = (1.3 ± 1 .2)% (to the BAT flux), /V/, = 
4.8! 2 1 5 9 x 10 22 cm" 2 (49+ 3 2 5 4 %), and TV 2 = 1.9/° 0 4 9 x 10 22 cm' 2 
(51/. 2 3 5%). They are all consistent with the case of = 100 r g 
within the statistical errors at the 90% confidence level. 

3.2.3. Narrowband Spectral Analysis of Iron-Ka Line Features 

To confirm the results on the reflection components obtained 
from the analysis of the broadband spectra, we also analyze the 
narrowband spectra of the XIS-FIs and XIS-BI in the 3-9 keV 
and 3-8 keV bands, respectively, focusing on the fluorescence 
iron-Ka line. Figure 6 plots the spectra of the XIS-FIs. We model 
the continuum by the same model from which iron-K emission 
lines are excluded by replacing the pexmon model with pexrav. 
A narrow Gaussian and a diskline component (Fabian et al. 
1989) are added, representing iron-K line emission from the 
torus and accretion disk, respectively. Although iron-K lines 
from radio-loud AGNs are usually narrow (e.g., Sambruna et al. 
2009; Larsson et al. 2008) and hence are likely produced by the 
torus, here we aim to constrain any possible contribution from 
the accretion disk directly from our data. Both line energies are 
fixed at 6.4 keV in the rest frame. In the diskline model, the 
inner and outer radii are fixed at r m = 1 00 r„ and r out = 10 5 r g , 
respectively, with an emissivity law of r~ 3 . The parameters 
of the continuum are fixed at those in Table 2 except for the 
normalization. 

The fitting results are summarized in Table 3. Data to 
continuum model ratio is shown in the lower panel of Figure 6. 
We find that the observed equivalent widths of the diskline 
component are < 53 eV (3C 206) and < 68 eV (PKS 0707-35), 
and those of the narrow line are < 30 eV (3C 206) and 58±26 eV 
(PKS 0707—35), respectively. To discuss the physical origin 
of these lines, we need to correct the observed equivalent 
widths for the contribution of the jet in 3C 206. The jet flux 
accounts for ~26% of the total flux at 6.4 keV in the rest 
frame. We thus obtain the equivalent widths of the broad line 
EW d"k < 71 eV and the narrow-line EW^ SS < 40 eV for 3C 
206. Furthermore, for the line originating from the torus, it is 
necessary to take into account the time variability of the direct 
component. The equivalent widths of the narrow iron-K line 
with respect to the five-year averaged BAT flux are found to be 
EW“ r2 ss < 71 eV and 72±36eV for 3C 206 and PKS 0707—35, 
respectively, which can be directly compared with theoretical 
predictions described in Sections 3.2.1 and 3.2.2. These results 
are consistent with those obtained from the broadband spectra 
in both targets. 

When we assume the inner radius of the accretion disk to be 
10 r g , the equivalent widths corrected for the jet contribution and 
time variability become EW“ a r2 ss = 42 ± 38 eV and 79 ± 33 eV 
and HW™' 4 < 67 eV and 1 15 ± 64 eV for 3C 206 and 
PKS 0707—35, respectively. These values are still consistent 
with those derived for the case of r m = 1 00 r g . However, the 
equivalent width of the diskline component in 3C 206 is smaller 
than the predicted value of 100-130 eV by George & Fabian 
(1991; see Section 3.2.1). Recalling the fact the model with 


3C 206 




Figure 6. (a) Folded XIS-FI spectrum of 3C 206 in the 3-9 keV band (black 
crosses). In the spectral fitting, we include the XIS-BI spectrum in the 3-8 keV 
band, which is not plotted here for clarity. The top solid lines show the best-fit 
continuum model in Table 2. The best-fit narrow and broad iron-Ka lines are 
represented with the red solid and blue dotted lines, modeled by a Gaussian 
and a diskline profile, respectively. The ratio between the data and continuum 
model is plotted in the second panel. Note that the feature around 7.5 keV is 
instrumental due to a nickel-Ka line in the background spectrum (Tawa et al. 
2008). (b) The same for PKS 0707-35. 

(A color version of this figure is available in the online journal.) 


r ln = 1 00 r g better reproduces the broadband spectrum, we 
infer that the accretion disk of 3C 206 is likely to be truncated 
at a radius much larger than 10 r g if the disk produces only cold 
reflection. 


3.3. Constraining the Torus Structure 

Assuming that the narrow iron-Ka line entirely originates 
from the torus, we can constrain the solid angle and/or hydrogen 
column density of the torus from the line equivalent width. To 
quantify the relation between the predicted equivalent width of 
the iron-K line and torus geometry, we employ the numerical 
model by Ikeda et al. (2009), where the reprocessed emission 
from the torus irradiated by the central source is calculated via 
Monte Carlo simulations. 

Ikeda et al. (2009) consider a three-dimensional torus con- 
figuration, characterized by the half-opening angle 0 oa , the in- 
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e oa ( de 9) 

Figure 7. Expected equivalent width of the iron-Ko: line to the total continuum 
given as a function of the opening angle # oa , calculated from the Monte Carlo 
simulation by Ikeda et al. (2009). The results are given for two equatorial column 
densities, N^ q = 10 23 (filled marks) and 10 24 cm' 2 (opened marks). The circles 
(black), squares (red), and triangles (blue) represent the results at the inclination 
angle #j nc = 8°, 28°, and 48°, respectively. In the range between $ oa = 70° and 
90°, the equivalent widths are extrapolated by assuming that the line intensity 
is proportional to the volume of the torus. The dashed lines (magenta) show the 
90% upper and lower limits of the equivalent width obtained from the spectra 
of 3C 206 when rj n = 100 r g is assumed for the disk-reflection component (see 
the text). 

(A color version of this figure is available in the online journal.) 


Table 3 

Best-fit Parameters of Suzaku Narrowband Spectra 




3C 206 

PKS 0707-35 

(1) 

*§auss ( keV ) 

6.4 a 

6.4 a 

(2) 

cr (eV) 

l a 

l a 

(3) 

<keV> 

6.4 a 

6.4 a 

(4) 

Lin (Eg) 

100 a 

100 a 

(5) 

Afeauss (photon cm' 2 s” 1 ) 

(<4.3 x 10“ 6 ) 

(4.7 ±2.1) x 10“ 6 

(6) 

EWg> a s uss (eV) 

(<30) 

58 ±26 

(7) 

EW£ 2 . (eV) 

(<71) 

72 ±36 

(8) 

A/disk (photon cm -2 s -1 ) 

{<1.1 x 10' 6 ) 

(<5.4 x 10" 6 ) 

(9) 

EW^ k (eV) 

(<53) 

(<68) 

(10) 

EW“ r ‘ (eV) 

(<71) 

(<68) 


X 2 /dof 

312.50/304 

389.33/366 


Notes. The results of narrowband spectral fitting (3-9 keV for XIS-FIs and 
3-8 keV for XIS-BI). The errors correspond to 90% confidence level for a single 
parameter. If we derive only the upper limit, it is shown in the parenthesis. 

(1) The center energy of the narrow iron-Ko line in the rest frame. 

(2) The line width for the narrow iron-Ko line. 

(3) The center energy of the broad iron -Ko line in the rest frame. 

(4) The inner radius of the accretion disk. 

(5) The total photon flux of the narrow iron-Ko line. 

(6) The observed equivalent width of the narrow iron-Ko line with respect to 
the Suzaku continuum flux. 

(7) The corrected equivalent width of the narrow iron-Ko' line with respect to 
the time-averaged continuum flux measured with .S’ re//// I1AJ" after subtracting 
the jet component (in the case of 3C 206). These equivalent widths are written 
as EW“ a r2 ss in the text (see Section 3.2.3). 

(8) The total photon flux of the broad iron-Ko line. 

(9) The observed equivalent width of the broad iron-Ko line with respect to the 
Suzaku continuum flux. 

(10) The corrected equivalent width of the broad iron-Ko line with respect to 
the Suzaku continuum flux after subtracting the jet component (3C 206 only). 
This equivalent width is written as EW/° r k in the text (see Section 3.2.3). 

a The parameters are fixed. 


o 



Figure 8. Same as Figure 7, but calculated for type-2 AGNs where the inclination 
is larger than the opening angle of the torus. /V^ q = 1 O 23 cm' 2 and A inc = 88° are 
assumed. The dashed lines (magenta) show the 90% upper and lower limits of the 
equivalent width obtained from the spectra of PKS 0707—35 when r m = 1 00 r g 
is assumed for the disk-reflection component (see the text). 

(A color version of this figure is available in the online journal.) 

clination angle of the line-of-sight /9 mc , the hydrogen column 
density in the equatorial plane /V^ q , and the ratio between the 
inner radius r m and outer radius r out , which is fixed at 0.01 
(see Figure 2 of Ikeda et al. 2009). A cutoff power law with a 
cutoff energy of 360 keV is assumed for the spectrum of the 
central source; this choice little affects the predicted iron-K line 
intensity when compared with the case of 200 keV. For a pho- 
ton index of T — 1.8 and a column density of /V,^ 4 = 10 23 or 
10 24 cm -2 , we calculate the iron-K line equivalent widths for 
various combinations of 0; nc and 0 oa , using the spectral database 
available as a “table model” in XSPEC. The results applicable 
for type-1 AGNs (i.e., 0i nc < 6 0 a ) and for type-2 AGNs (i.e., 
$inc > #oa) are displayed in Figures 7 and 8, respectively. The 
inclination is fixed at $i nc = 8°, 28°, or 48° for type-1 AGNs 
and at 0j nc = 88° for type-2 AGNs. 8 As representative values, 
we assume A^ q = 10 23 or 10 24 cm -2 in Figure 7, and A', 1 ) 4 = 
10 23 cm -2 in Figure 8, the same order of the line-of-sight col- 
umn density of PKS 0707—35. Although the assumed value of 
/Vy 4 more sensitively affects the observed equivalent width in 
type-2 AGNs than in type-1 AGNs because of absorption of the 
continuum, the difference in the attenuation fraction of the direct 
component at 6.4 keV between the line-of-sight column density 
of 3 x 10 22 cm~ 2 (as observed) and 10 23 cm 2 is found to be 
~10% level, which we ignore in the following discussion for 
PKS 0707—35. Since the Ikeda et al. (2009) model is available 
only for 0 m < 70°, we extrapolate the result of 0 oa = 70° toward 
0 m = 90° in Figure 7, assuming that the line intensity is pro- 
portional to the volume of the torus. The figures show the trend 
that the equivalent width decreases with increasing half-opening 
angle except for the case of /V^ q = 10 24 cnT 2 in Figure 7, which 
is consistent with Figure 13 of Ikeda et al. (2009). 

In Figures 7 and 8, we also show the range of the “corrected” 
equivalent width of the narrow iron-Ka of 3C 206 and PKS 
0707—35, respectively, with dashed lines. We find that the torus 
of 3C 206 either has a relatively low column density in the 
equatorial plane /V^ q < 10 23 cm' 2 , or has a large opening angle 
(0 O a > 80°)if/V^ q = 10 24 cnT 2 . Even when we assume 0.5 times 
smaller elemental abundances for iron and nickel than the solar 


The choice of a lower 0 mc value for type-2 AGN little changes the result. 
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values, similar constraints are obtained; iV^ q < 10 23 cm 2 , or 
$ oa > 70° for (V^ 1 = 10 24 cm' 2 . The result indicates that the 
torus of 3C 206 is far from Compton-thick or covers only a 
small solid angle. The line-of-sight column densities of PKS 
0707—35, Nn = (1-8) x 10 22 cm -2 , suggest a small A^ q for 
this object, and we confirm that the iron-Ka equivalent width is 
consistent with the prediction from (V^ q ~ 10 23 cm' 2 as shown 
in Figure 8. 

4. SUMMARY AND DISCUSSION 

4.1. Results Summary 

We have obtained the best quality broadband X-ray spectra 
of the two most luminous radio-loud quasars 3C 206 (type-1) 
and PKS 0707—35 (type-2) with Suzaku and Swift/ BAT, which 
have hard X-ray luminosities of 10 45 5 erg s' 1 and 10 44 9 erg s' 1 
(14-195 keV), respectively. The spectra are well described 
with a physically motivated model consisting of the direct 
component (with partial absorptions for PKS 0707—35), its 
reflection component (including fluorescence lines) from the 
accretion disk, that from the torus, jet emission (3C 206), 
and scattered component (PKS 0707—35). Here, we properly 
take into account the effect of time variability of the direct 
component between the Suzaku and (5 yr averaged) Swift /BAT 
observations. Our main findings are that (1) there is a significant 
contribution from the jet in the X-ray spectrum only for 3C 206 
and that (2) the reflection strength from the torus is estimated 
to be quite small in both objects compared with typical Seyfert 
galaxies, which is confirmed by the detailed analysis of the 
iron-K emission line. 

4.2. Jet Contribution to the X-Ray Spectra 

We have shown that the X-ray spectrum of the “type-1” 
radio-loud quasar 3C 206 is likely to contain significant emis- 
sion from the jet, while its contribution is negligible for the 
“type-2” radio quasar PKS 0707—35. In the unified scheme of 
AGNs, radio-loud quasars and radio galaxies are intrinsically 
the same as blazars except for the inclination angle with respect 
to the jet axis, which is the smallest in blazars and increases 
from type-1 to type-2 AGNs. We show that it is possible to 
make a rough estimate of the contribution from the jet in the 
X-ray band using the radio core luminosity, by referring to the 
“blazar sequence” with correction of beaming factor based on 
reasonable assumptions. Our results imply that non-thermal jet 
emission may make a significant contribution to the observed 
X-ray spectrum of type-1 radio-loud AGNs in general, and this 
component should be taken into account when spectral modeling 
is done. 

The broadband spectral fit suggests that the fraction of the jet 
component in the X-ray spectrum of 3C 206 is 26% (2-10 keV) 
and 62% (14-195 keV) of the total flux, which is consistent with 
the estimate from the radio luminosity within a factor of two 
If this is the case, the very large hard X-ray luminosity of this 
source is partially due to the jet emission. The same trend that 
the relative jet power increases with energy is also reported from 
3C 273 (Grandi & Palumbo 2004). The true luminosity arising 
from the accretion disk is estimated to be 2.2 x 1 0 44 erg cm' 2 s' 1 
(2-10 keV) and 1.2 x 10 45 erg cm' 2 s' 1 (14-195 keV), which 
is very similar to that of PKS 0707—35. 

4.3. Reflection Components and Torus Structure 

Compton reflection components accompanied by fluores- 
cence lines from matter around the nucleus give us valuable 


information on the geometry of the surroundings of an AGN. 
Unlike in most of previous works, we separately consider two 
reflection components, one from the accretion disk and the other 
from the torus in order to achieve a physically consistent pic- 
ture. Relativistic blurring is considered in the disk reflection, 
whose parameters are linked to the reflection strength in a phys- 
ically consistent manner, based on the calculation by George 
& Fabian (1991). Since the torus is located far away from the 
central SMBH and CCD spectroscopy cannot resolve the line 
width, we ignore any blurring for the latter. It is hard to dis- 
tinguish the two reflection components from the spectra alone, 
which would be strongly coupled with each other in the fitting. 
Thus, we assume the reflection strength from the accretion disk 
to be R d isk = 0.1 (3C 206) and 0.2 (PKS 0707-35), and then 
derive the reflection strengths from the torus as a free parameter. 
We obtain R toms = 0.29 ±0.18 (3C 206) and 0.41 ± 0. 18 (PKS 
0707—35) from the broadband spectral fitting. These results are 
consistent with the equivalent widths of the iron-Ka line derived 
from narrowband (3-9 keV) spectral analysis. 

Since any possible broad iron-K line features look very weak 
in both objects, it is difficult to directly constrain the innermost 
radius of the accretion disk from the line profile. We find that the 
broadband fitting of the 3C 206 spectra becomes significantly 
worse if we assume r m = 1 0 r g , which predicts a much larger 
reflection strength from the disk, compared with the case of 
ri n = 1 00 r„. In addition, the iron-K line analysis in the narrow 
band yields an upper limit of the equivalent width of a broad 
iron-K line that is smaller than the prediction by George & 
Fabian (1991) in the case of r ln = 10 r g . These results imply 
that the accretion disk of 3C 206 is likely to be truncated at 
a radius much larger than r m = 1 0 r g as far as the disk is not 
significantly ionized. In fact, if we make the disk inner radius a 
free parameter, by assuming the relation between r m and R based 
on the equivalent width calculated by George & Fabian (1991) 
and its corresponding reflection strength in the pexmon model 
(Nandra et al. 2007), then we obtain the best-fit inner radius of 
r m = 50 49 o r g . Pictures of truncated disks are obtained from 
other type-1 radio-loud AGNs, 4C 50.55 (r ln > 340 r g ; Tazaki 
et al. 2010), 3C 111 (r in ~ 20-100 r g ; Tombesi et al. 2011), 
3C 390.3 (r in > 20 r g ; Sambruna et al. 2009), and 4C+74.26 
(/- m > 44 r„ ; Larsson et al. 2008). By contrast, the disks of 3C 
120 and 3C 382 seem to be extended down to r ln = (9 ± 1) r„ 
(Kataoka et al. 2007) and r; n = 1 2 ± 2 r„ (Sambruna et al. 2011), 
respectively. It is not clear if there is any systematic trend of disk 
truncation in radio-loud AGNs. We note that the presence of a 
highly ionized inner disk extending to r m = 6 r g , or smaller if 
the black hole is rapidly spinning, cannot be ruled out, however, 
which might explain the soft excess seen in 3C 206. 

Our robust finding is that the reflection component from the 
torus containing a narrow iron-K line is weak in both 3C 206 
and PKS 0707—35 (R torus < 0-6) compared with the typical 
reflection strength observed in the X-ray spectra of Seyfert 
galaxies (R ~ 1; Dadina 2008). Utilizing the torus model by 
Ikeda et al. (2009), we have quantified the relation between 
the half-opening angle of a torus (0 oa ) and the equivalent 
width of an iron-K line. The equivalent width with respect 
to the time-averaged, non-jet continuum flux is < 71 eV for 
3C 206. This result constrains that the column density in the 
equatorial plane A^j q < 10 23 cm' 2 , or 0 oa > 80° if A^ q = 
10 24 cm' 2 is assumed, which corresponds to a very small 
solid angle of Q torus /4^ < 0.17. Also, a small value of lV H q , 
~10 23 cm' 2 , is inferred for PKS 0707—35 from the line-of- 
sight absorption column density, which can explain the iron-K 
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line equivalent width. Thus, these luminous radio-loud AGNs 
have only poorly developed tori. There is a possibility that the 
tori are geometrically thick but highly clumpy with a very small 
filling factor. 

It is interesting to compare this result with those obtained 
from less luminous radio-loud AGNs. Eguchi et al. (2011) 
analyze Suzaku and Swift /HAT spectra of the radio galaxy 
NGC 612, which has Z. 2-10 = 3 x 10 43 erg s' 1 , and obtain 
the equivalent width of a narrow iron-K line of 280 ± 60 eV. 
They also find that the torus model by Ikeda et al. (2009) 
well represents the spectra with the best-fit parameters of 
A^ q ~ 1.1 x 10 24 cm' 2 and 6 m ~ 60°-70°. The opening 
angle of NGC 612 is smaller than the estimate for 3C 206, 
0 O a > 80°, obtained when a similar value of A^j q = 10 24 cm' 2 
is assumed. Two narrow-line radio galaxies studied by Tazaki 
et al. (2011), 3C 403 (L 2 -io = 2 x 10 43 erg s' 1 ) and IC 
5063 (L 2-10 = 6 x 10 42 erg s' 1 ), show narrow iron-K line 
with equivalent widths of 460 ± 100 eV and 150 ± 20 eV, 
respectively. Thus, the iron-K equivalent widths obtained from 
type-2 radio-loud AGNs with moderate luminosities are all 
larger than that of most luminous one PKS 0707—35 (L 2-10 = 
2.1x 1 0 44 erg s' 1 ). This can be partially explained by their larger 
line-of-sight column densities that attenuate the continuum 
fluxes at 6.4 keV: JV H = 1.1 x 10 24 cm' 2 (NGC 612), 
6.1 x 10 23 cm' 2 (3C 403), and 2.5 x 10 23 cm' 2 (IC 5063). 

The comparison suggests that the torus solid angle and/or 
column density are smaller in more luminous radio-loud AGNs. 
It is known that the fraction of absorbed AGN decreases with 
X-ray luminosity above L 2-10 ~ 10 42 erg s' 1 from studies 
using mostly radio-quiet AGN samples (e.g., Ueda et al. 2003; 
Hasinger 2008; Burlon et al. 2011). This fact indicates that 
the solid angle of the torus decreases with luminosity in radio- 
quiet AGNs. A popular explanation is the receding torus model 
(Lawrence 1991), where the strong radiation from the central 
engine affects the structure of the torus; both the solid angle and 
column density could be reduced in luminous AGNs. Our results 
are consistent with the idea that the same picture holds for radio- 
loud AGNs. This implies that the presence of powerful jets is not 
related to the structure of the torus, which is responsible for gas 
feeding into the disk as the boundary condition. Rather, some 
(not yet established) physics taking place in more inner parts 
of the disk must be involved to explain the difference between 
radio-loud and radio-quiet AGNs. 
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